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ABSTRACT 

Supermassive black hole binaries (BHBs) produced in galaxy mergers recoil at the time of their coalescence 
due to the emission of gravitational waves (GWs). We simulate the response of a thin, two-dimensional disk 
of collisionless particles, initially on circular orbits around a 10'' M0 BHB, to kicks that are either parallel 
or perpendicular to the initial orbital plane. Typical kick velocities (vkick) can exceed the sound speed in a 
circumbinary gas disk. While the inner disk is strongly bound to the recoiling binary, the outer disk is only 
weakly bound or unbound. This leads to differential motions in the disturbed disk that increase with radius 
and can become supersonic at > 700 Schwarzschild radii for yyck = 500 km s"', implying that shocks form 
beyond this radius. We indeed find that kicks in the disk plane lead to immediate strong density enhancements 
(within weeks) in a tightly wound spiral caustic, propagating outward at the speed ^ vyck- Concentric density 
enhancements are also observed for kicks perpendicular to the disk, but are weaker and develop into caustics 
only after a long delay (>one year). Unless both BH spins are low or precisely aligned with the orbital angular 
momentum, a significant fraction (> several %) of kicks are sufficiently large and well aligned with the orbital 
plane for strong shocks to be produced. The shocks could result in an afterglow whose characteristic photon 
energy increases with time, from the UV lOeV) to the soft X-ray lOOeV) range, between one month and 
one year after the merger. This could help identify EM counterparts to GW sources discovered by LISA. 
Subject headings: black hole physics - galaxies: nuclei - gravitational waves 



1. INTRODUCTION 

The recent break-through in numerical relativity has al- 
lowed a direct computation of the linear momentum flux 
produced during the coalescence of a BH binary (Baker et 
al. 2006, 2007; Campanelli et al. 2007a,b; Gonzalez et al. 
2007a,b; Herrmann et al. 2007a,c; Koppitz et al. 2007). The 
resulting final recoil depends on the masses, orbital parame- 
ters, and spins of the BHs, and in special configurations (i.e. 
with spins anti-aligned with each other), it can reach veloci- 
ties as high as « 4,000 km s"' (Gonzalez 2007b). Much of 
the recent literature focused on the astrophysical implications 
of high-velocity kicks, which may displace or remove super- 
massive BHs from galactic centers (e.g. Merritt et al. 2004; 
Madau & Quataert 2004; Loeb 2007), and inhibit the growth 
of BHs at high redshifts (z > 6), where the escape velocities 
from low-mass galactic halos is small (< 100 km s"' ; Haiman 
2004, Yoo & Miralda-Escude 2004, Shl^piro 2005; Volonteri 
& Perna 2005). 

Another exciting implication of kicks is that they may help 
produce an electromagnetic (EM) counterpart of gravitational 
wave sources detected by the future Laser Interfewmetric 
Space Antenna {LISA) satellite. The discovery of such a coun- 
terpart would constitute a milestone for fundamental physics 
and astrophysics (e.g. Kocsis et al. 2007b). If the BHB is 
surrounded by a circumbinary gas disk, the disk will respond 
promptly (on the local orbital timescale) to such a kick. If this 
results in warps or shocks, the disturbed disk could produce a 
transient EM signature (Milosavljevic & Phinney 2005). The 
sky localization uncertainty from the LISA instrument several 
weeks prior to merger is typically a few square degrees (Koc- 
sis et al. 2007a; Lang & Hughes 2007), and the kick can begin 
building up at the end of the inspiral phase, before the final 



coalescence (Schnittman et al. 2007). This will make it pos- 
sible, in many cases, to monitor a few square-degree area on 
the sky prior to, during, as well as immediately following the 
coalescence of BHs in the mass range ~ ( 10^-10^) M0/(l+z) 
at redshifts out to z ~ 3, and search for a prompt transient sig- 
nature associated with the kick (Kocsis et al. 2007b). 

In the context of producing a prompt EM counterpart, we 
expect that the direction, in addition to the magnitude of the 
kick, will be important. Naively, one expects that a kick 
within the plane of any circumbinary disk will be more likely 
to cause density enhancements and "light up" the disk than 
a kick perpendicular to it. The kick direction can be clearly 
important on larger scales, as well, and affect phenomena that 
would occur on long timescales {':$> years) after the merger 
For example, the angle relative to a large-scale galactic disk 
can determine whether a kicked BHB ends up outside the 
galaxy or not, and perhaps also whether shocks are produced 
when a BHB plunges into such a large-scale disk. While 
these issues merit investigation, in this Letter we focus on the 
smaller scales and shorter timescales (< months) that are rel- 
evant to prompt LISA counterparts. 

Schnittman & Buonanno (2007; hereafter SB07) used the 
"effective one-body approach" and derived a scaling formula 
that yields the recoil velocity vector for arbitrary mass ratios 
and spin vectors. While their results have not been tested 
for generic spins, they agree well (to within 20-30%) with 
numerical results in those special configurations where they 
were tested (including numerical calculations for configura- 
tions with the spins parallel or at intermediate angles with 
respect to the orbital plane; see Campanelli et al. 2007a,b, 
Gonzalez et al. 2007a; Herrmann et al. 2007b, Tichy & Mar- 
ronetti 2007). 
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In this Letter, we investigate the response of a circumbi- 
nary disk to the kick, by following the perturbed, Keplerian 
orbits of collisionless massless test particles around the recoil- 
ing BHB. The purpose of this exercise is to demonstrate that 
prompt shocks or strong density enhancements are likely to 
arise when the kick is aligned with the plane of the circumbi- 
nary disk, whereas they may be less likely for highly inclined 
kicks. We will then use the formula of SB07 for the depen- 
dence of the kick speed and direction on the mass ratio and 
spins of the BHs, to argue that a significant fraction (> few 
%) of kicks may be both sufficiently large and well aligned 
with the orbital plane for strong shocks to be produced within 
a few weeks after coalescence. 

2. CIRCUMBINARY DISKS 

We begin with the assumption that the kicked BHB is 
surrounded by a rotationally supported, geometrically thin 
gaseous circumbinary disk. The idealized case of a smooth 
axisymmetric accretion disk, aligned with the binary's orbital 
plane (Bardeen & Peterson 1975), with a simple vertical struc- 
ture can be described by its density and temperature profiles 
(p(r), T(r)), as well as the scale height //(r). This minimal in- 
formation is needed (i) to compute the angle extended by the 
disk as viewed from the center, and (ii) to examine whether 
the kick is supersonic (v^ck > c^; a necessary condition for 
shocks to occur). In addition, we need to know the overall 
disk size (necessary to assess what fraction of the disk gas 
remains bound to the recoiling BHB), and whether the disk 
mass is comparable to the BHB mass (necessary to assess 
whether the BHB is slowed down by the disk). 

Unfortunately, fully self-consistent and stable accretion 
disk models around supermassive BHs are both difficult to 
produce, and require many ad-hoc assumptions (see, e.g., 
the review by Blaes 2007 for a list). While variants of 
the parametric so-called a-disk models (Shakura & Sun- 
yaev 1973) have been successfully calibrated against ob- 
servations of stellar-mass objects, observations of active 
galactic nuclei (AGN) have not led to similarly robust con- 
straints on disk models around supermassive BHs. Never- 
theless, several authors have di scussed the p ossible be havior 
of a circumbinary g-disk (e.g.lArmitage & Nataraian 120021 : 
iMilosavlievic & PhinnevI 120051: iDotti et all 12006) " A key 
quantity is the critical orbital semi-major axis Ucm at which 
the time-scale for the decay of the orbit due to gravitational 
radiation becomes shorter than the time-scale for viscous 
evolution (Begelman et al. 1980). Once the binary separa- 
tion shrinks below this value at fciit, the gas outside will not 
evolve significantly before the BHs coalesce. Most of the 
"kick" momentum is accumulated during the final stages of 
coale scence, on a rapid time-scale dSchnittman & Buonanno I 
l2007h . hence a snapshot of the disk at fcrit outside ^ acrft 
should still describe the disk around the kicked binary. Prior 
to fcrit, the torques from the binary may create a central cav- 
ity nearly devoid of gas within the radius rcavitv ^ 2acrit (for 
a nearly equal-mass binar y, e.g.^A rtymowicz & Lubowll994t 
iMacFadven & Milosavliev ic 2006), or a narrower gap around 
the orbit of the lower-mass BH in the case of unequal masses 
q=My /M2 <C 1 (e.g. Armitage & Natarajan 2002). In the lat- 
ter case, the lower-mass hole will "usher" the gas inward as its 
orbit decays, so by the time of coalescence, this gas may ac- 
crete onto the more massive BH, again leaving a nearly empty 
central cavity of radius Ocrit- 

In the inner regions of a disk around a supermassive BH, the 
dominant vertical support is expected to be radiation pressure 



(as opposed to gas pressure), and ele ctron scattering domi- 
nates the opacity over absorption (e.g. Shapiro & TeukolskyJ 
I1983L p. 440). Milosavlievic & Phinney (2 005 ) presented an 
explicit model for a circumbinary disk with a central cavity, 
for the case of a nearly equal-mass binary at fcrit- They find, 
under the above conditions for the inner disk (rcavity = 2acrit), 

rcavtty=117rsa:?•^^ry_,/m)«•2X■«'[4?/(l+?)']°•'^ 

where Mg denotes the total BHB mass M\ + M2 in units of 
IO'^Mq, rs = 2GM/c^ is the Schwarzschild radius, O.la-i is 
the effective a-parameter relating the kinematic viscosity to 
the gas (not total) pressure, m is the mass accretion rate at 
'"cavity in units of the Eddington rate Mndd (i-C- the accretion 
rate that would produce the Eddington luminosity with the 
radiative efficiency 0.1 77-1), and q = M\/M2 < 1 is the mass 
ratio (note that we neglect here the dependence on additional 
dimensionless coefficients of order unity). 
The thickness of the disk at its inner edge rcavity is given by 

^ = 0.46M,-°'2[4^/(l +^)2]-l-84^0.76(„./^_^)2.43 

and the mid-plane gas temperature is 
r= 1.7 X 10X''-2'[4?/(l+?)']"''"'«°f (m/^-i)°''K. (2) 
The ratio of radiation and gas pressure is given by 

^gas 

Note that the disk is marginally thin at rcavity The expected 
radial profile of such a disk outside rcavity was discussed by 
Goodman & Tan (20 0^. In th e stable case considered by 
IMilosavlievic & PhinnevI ( 1200 5), when the viscosity is pro- 
portional to the gas (rather than the total) pressure, the mid- 
plane temperature and density, and surface density vary with 
radius as T oc r"^/'", p oc r~^l^, and E cx r~'^l^ . The ratio of 
gas/total pressure varies as /3 = Pgas/Aot oc r^'/'" (in the limit 
of /3 ^ 1). Inside the transition radius where /? « 1, 

r,, « 1.9 X 10X^"af/,^i(m/'?-i)''/'''-s, (4) 

the scale height h is constant, or /i/r oc 1/r (here atot-i relates 
the kinematic viscosity to the total pressure). This further as- 
sumes that the disk is optically thick. Beyond rn, gas pres- 
sure becomes dominant for vertical support, and the scale- 
height begins to rise as /z cx r^Vzo Somewhat farther out, 
beyond ^ 10,000rs (for Me ^ 1), the disk will become self- 
gravitating and subject to Toomre-instability. There are both 
theoretical and observational reasons to suspect that disks in 
real AGN extend to larger radii (Goodman 2003; Blaes 2007). 
However, note that the orbital time-scale at this radius is «2 
years, so that in the context of prompt LISA counterparts, we 
are not interested in radii outside this annulus. 

To summarize, based on the above, we adopt the follow- 
ing simplified picture for the disk around a fiducial M\ +M2 = 
10^ Mq binary. The disk has an inner edge at lOOrs (inside 
which it is empty) and an outer edge at 10,000rs (outside 
which there may still be gas, but it evolves slowly and we 
do not follow it). The scale-height and temperature at the in- 
ner edge is /i/r = 0.46 and T = 1.7 x 10''K, respectively. The 
scale-height remains constant with radius out to 2, OOOrs, be- 
yond which it increases nearly linearly (h cx r^'/^*^). The tem- 
perature varies with radius as T cx r"^/'". 



Fig. 1 . — The top view of a disk around an M[ +M2 = lO' Mq BH binary, 
which recoiled within the disk plane at the velocity of v^g^^i^ = 500 km s"' , ori- 
ented vertically upward in the diagram. The disk is initially assumed to have 
an inner edge at r|„ = 1 OOrg , and is shown here out to a radius of = 5 , OOOrj . 
at a time f = 90 days after the kick. A similar tightly wound spiral caustic pat- 
terns develops already in less than a month. The dark/light shades con'espond 
to regions of high/low density, the low-density regions being similar to the 
initial surface density, and the high-density regions about 10 times overdense. 

The important features of such a disk (as well as other pro- 
posed variants of a-disks) are the following: (i) orbital mo- 
tions in thin disks are supersonic (Pringle 1981), so that the 
gas is susceptible to shocks if disturbed; (ii) at the relevant 
radii outside lOOrs the viscous time-scale is long, so that the 
orbits are near Keplerian; (iii) gas near the inner edge of disk 
is tighdy bound to the kicked BHB (vorbit ~ 3 x lO"* km/s), but 
the outer edge (vorbit ^ 3 x lO-' km/s) can be marginally bound, 
or even unbound, for large kicks (the approximate condition 
for being bound is Voi-bit > 2.4vkick), and (iv) the total disk mass 
within lOjOOOrs is much less than the BHB mass, which jus- 
tifies ignoring the inertia of the gas bound to the BHB. 

3. THE RESPONSE OF THE DISK TO KICKS 

For a quantitative assessment of the disk's response to the 
kick, we employ the following approximation: the disk par- 
ticles are assumed to be massless, collisionless, and initially 
on co-planar, circular orbits. The kick simply adds the veloc- 
ity Vkick to the instantaneous orbital velocity of each particle 
(in the inertial frame centered on the BHB). We used = 10^ 
particles, distributed randomly and uniformly along the two- 
dimensional surface of the disk. The kick velocity was varied 
between 500 km s"' < v^ck < 4,000 km s"', and directed ei- 
ther perpendicular or parallel to the initial disk plane. 

Figure [T] shows, as an example, a face-on view of the sur- 
face density of the disk 90 days after a kick with vyck = 
500 km s"' in the plane of the disk (/kick = 0°). The sharp, 
tightly wound spiral features clearly seen in the figure trace 
the locus of points where particles cross each other, corre- 
sponding formally to a density caustic. The spiral caustic first 
forms at ^ 30days, and then propagates outward at a speed 
of « 500 km s"'. This behavior can be roughly understood as 
follows: at a given radius rc, the caustic forms at the time 
when the radial epicyclic motions from two neighboring an- 
nuli, separated by the epicyclic amplitude ^ (vkick/i'orbit)'"c, 
overlap (this predicts ^ fci^'kick)- We also found that as Vkkk 
is increased, the spiral caustics spread farther out. Eventu- 
ally, for kicks strong enough so that a significant fraction of 
the disk is unbound, the caustics loose their coherent spiral 
patterns, and develop complex two-dimensional shapes. 

Figure |2] shows a side view of the 3D particle density one 
week after a kick with the velocity Vkkk = 500 km s"' perpen- 
dicular to the disk (/kick = 90°). The density profile in this 
case remains azimuthally symmetric, but still develops con- 
centric rings of density fluctuations. The important difference 
from the parallel kick case is that the density enhancements 



Fig. 2. — The aerial view of a disk as in Figure[T] except here the kick is 
oriented perpendicularly to the disk, and the snapshot is taken at ? = 1 week. 
For visual clarity, the graphic contrast was increased relative to Fig.[T] and the 
diagram stretched by a factor of 10 along the axis perpendicular to the disk 
plane. The density distribution is azimuthally symmetric, and while there 
are mild concentric density fluctuations, strong enhancements (i.e. caustics) 
develop only after a delay of ss one year 

are much weaker (at the ten percent level). By examining 
the time-evolving radial cross-sections, we have verified that 
sharp density enhancements, i.e. true caustics caused by the 
orbit-crossing of particles, first appear only after one year, 
and involve a smaller fraction of the disk particles. 

More generally, kicks can be expected to occur in a direc- 
tion intermediate between the two extremes shown in Figs. [T] 
and|2l Following a kick at an arbitrary inclination angle, each 
particle still remains on a Kepler orbit. Each new orbit will 
follow an ellipse in a plane that is tilted by an angle of order 
between ivkick/i^'orbit about the axis connecting the instanta- 
neous position of the particle with the central BHB. Hence, 
after an orbital time, the particles of the disk will be smeared 
vertically, effectively thickening the original disk. The two- 
dimensional simulations we performed do not conclusive tell 
us a critical kick inclination angle for prompt density caus- 
tics to be produced - this would require three-dimensional 
simulations, following the orbits of particles in a 3D disk of 
finite thickness. However, since the prompt, strong density 
enhancements appear to develop in the plane of the kick, near 
the inner edge of the disk, one may conjecture that this critical 
inclination angle is given roughly by /yck < arctan(/i / r) (with 
h/r= 0.46 evaluated at rcavity)- 

In Figure [51 we show the distribution of kick angles, taken 
from SB07, assuming that both spins are randomly oriented. 
We require further that the component of vyck within the or- 
bital plane exceed 3 - 500 km s"' (« 3 - 5 x the sound speed; 
see justification in §|4] below). A strong kick within the or- 
bital plane is produced if the spins are large and parallel to the 
angular momentum, but anti-aligned with each other 

The main conclusion to draw from Figure[3]is that caustic- 
producing kicks are not rare - in the cases shown in the panels, 
they occur at least a few percent of the cases. This conclusion 
remains true unless (i) the spins of the BHs is significantly 
below the maximal Kerr value (a < 0.9), or (ii) the spins are 
significantly aligned with each other The latter may be ex- 
pected if th e spin of both BHs resul ts from accretion from the 
same disk ( Bogdanovic et alJl2007h . 

The criteria defined above for producing caustics may turn 
out very conservative. For example, prompt caustics could 
develop when particles from different layers of the disk cross 
along their orbits. A plausible weaker criterion for shocks is 
that the kick-induced tilts of the orbits do not thicken the disk 
beyond its original scale-height. Apart from trigonometric 
factors of order unity that describe its azimuthal dependence, 
the tilt angle is given by /tin ^ (v'kick/voibit)'kick; requiring only 
/tilt < arctan(/!/r), or /^ck < (vorbit/vkick)arctan(/j/r), density 
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Fig. 3. — The thick solid curves in both panels show the (arbitrarily nor- 
malized) differential probability distribution of kick inclinations relative to 
the disk plane, using the SB07 formula to compute vj^^k ■ Both panels assume 
an equal-mass binary = 1 ; we find results changed only by a few percent for 
q = 0.1). We assume random spin orientations. On the left panel we assume 
dimensionless spins of ai =ai = \, while on the right panel ai =02= 0.9. The 
dashed curves show the fraction of kicks at a given inclination angle that have 
a component in the disk exceeding a fixed threshold (vj^jj-i^ p > 500 km s"' in 
the left panel, and 300 km s"' in the right panel). The thin curves indicate 
the cumulative probability of sufficiently fast kicks with inclinations below a 
given value; the horizontal line shows the fraction of such kicks with inclina- 
tions < 24° (6 and 15% percent in the left and right panels). 

caustics could arise even for nearly perpendicular kicks. In 
follow-up work, we will employ three-dimensional simula- 
tions, to clarify the relevant caustic formation criterion. 

4. DISCUSSION 

The main result of this Letter is that strong density enhance- 
ments can form promptly after a supersonic kick in the plane 
of the circumbinary disk, within a few weeks of the coales- 
cence of a ~ 10^ M© BHB. Because the disk is cold, and 
caustics (Fig. [U are formed when particles first cross each 
other along their orbits, this implies that corresponding shocks 
could occur in a gas disk. For hydrodynamical shocks to occur 
within a finite-pressure gas, the relative motions Vc between 
the neighboring particles that produce the caustic must exceed 
the sound speed. At the outermost radius where the disk is 
marginally bound to the BHB, one expects Vc ~ Vkick VoMu 
relative motions will be slower further inside. The relative 
speed should roughly correspond to covering the epicyclic 
amplitude ^ (vkick/t^'oibit)'"c in the caustic-formation time t^ ^ 
'"f/vkick, yielding v,. - v^^|^/vorbit- For vyck = 500km s"', this 
predicts v'c ^ 25km s"'(r/1000rs)'''^; we have verified in our 
simulations that particles cross the caustics with speeds at 
about ^ 30% above this predicted value. Compared with the 
sound speed « 25km s"'(r/1000rs)~^''^'^, this suggests that 
the density waves produced by the kick in the gas beyond 
~ 700rs will indeed steepen into shocks. We also found that 



the inclination of the kick may be important in determining 
the strength and timing of such shocks - perpendicular kicks 
would only produce weaker density enhancements, at least 
until a delay of about a year A non-negligible fraction (> 
several %) of kicks could, however, be sufficiently large and 
well aligned with the orbital plane for shocks to be produced 
within a few weeks after coalescence. 

The nature of the emission resulting from the shocks or 
density enhancements will have to be addressed in future 
work, by computing the heating rate at the spiral shocks, 
and modeling the overall disk structure and vertical radi- 
ation transport. However, the disk's luminosity resulting 
from the kick, Lyck could be a non-negligible fraction of the 
BHB's Eddington luminosity, and therefore potentially ob- 
servable (Kocsis et al. 2007b). If Mshock = /shock(A^i +M2) 
is the mass of the shocked gas that is heated to tempera- 
tures corresponding to Vshock, and fshock is the time-scale on 
which the corresponding thermal energy is converted to pho- 
tons, then Lkick ~ (l/2)M,shockVshock/fshock, and Lkick/-LEdd ~ 
0.1(/shock/10-5)(yshock/500 kms-i)2(f,hock/lmonth)-i. We 
may also speculate on the spectral evolution of the "kick 
after-glow". Assuming Mshock oc Y^rdr cx r'^/"^ (with E cx 
r~^l^ and dr oc r'/^, the epicyclic amplitude), and Vshock ~ ^c, 
we find Lkick oc r^'*/'". This suggests that the luminosity may 
be dominated by the outermost shocked shells. The spectrum 
will then peak at the characteristic photon energy correspond- 
ing to A:rshock oc oc v^j^jijj oc r. The shocks could therefore re- 
sult in an afterglow, starting from rcavity/i^'kick ^ 30 days, first 
peaking in the UV band (^ lOeV), and then hardening to the 
soft X-ray (^ lOOeV) range after one year. The detection of 
such an afterglow would help identify EM counterparts to GW 
sources discovered by LISA. 

Our results need to be verified in three-dimensional simu- 
lations that resolve the orbits within a finite disk thickness. A 
realistic disk model, incorporating gas dynamics, is needed to 
study the correspondence between collisionless caustics and 
gaseous shocks. Finally, the SB07 formula needs to be con- 
firmed for generic spin configurations. Nevertheless, our re- 
sults do suggest that kicks due to gravitational waves may pro- 
duce a prompt EM signal. 

We thank A. Buonanno, M. Milosavljevic, K. Menou, and 
B. Kocsis for useful discussions, and M. Milosavljevic for 
suggesting the arguments based on epicycles. This work was 
supported by the Polanyi Program of the Hungarian National 
Office for Research and Technology (NKTH) and by NASA 
grant NNG04GI88G (to ZH). 
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